Dephasing of Cooper pairs and subgap electron transport in superconducting hybrids 
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We argue that electron-electron interactions fundamentally restrict the penetration length of 
Cooper pairs into a diffusive normal metal (N) from a superconductor (S). At low temperatures this 
Cooper pair dephasing length remains finite and does not diverge at T — >■ 0. We evaluate the 
subgap conductance of NS hybrids in the presence of electron-electron interactions and demonstrate 
that this new length can be directly extracted from conductance measurements in such structures. 



It is well known that a normal metal (N) attached 
to a superconductor (S) also acquires superconduct- 
ing propertiesii^. This superconducting proximity ef- 
fect is directly related to the phenomenon of Andreev 
reflections': At the NS interface Cooper pairs are con- 
verted into subgap quasiparticles (electrons) which can 
diffuse deep into the normal metal keeping information 
about a macroscopic phase of the superconducting con- 
densate. Such macroscopic quantum coherence of elec- 
trons in the normal metal gets destroyed by thermal fluc- 
tuations only provided the corresponding inverse elec- 
tron diffusion time (Thouless energy) becomes smaller 
than temperature T. As a result, superconducting co- 
herence extends into a normal metal at a typical length 
Lt ~ y/D/T (where D is the electron diffusion coeffi- 
cient) implying that the whole normal metal can demon- 
strate superconducting properties at sufficiently low T. 

This proximity induced superconductivity manifests it- 
self in a number of well known phenomena, such as Meiss- 
ner and Josephson effects in normal-superconducting 
hybrids'i'^ as well as dissipativc transport of subgap elec- 
trons across NS interfaces^. Provided the NS interface 
transmission is low its corresponding subgap (Andreev) 
conductance G remains rather small being proportional 
to the second order in the barrier transmission. On the 
other hand, G can be strongly enhanced at low ener- 
gies due to non-trivial interplay between disorder and 
quantum interference of electrons in the normal metal2rii 
which leads to the so-called zero-bias anomaly (ZBA) 
G oc 1/ VV and G oc I / \/T in the limit of low voltages 
and temperatures. 

In this paper we will demonstrate that in the low tem- 
perature limit both superconducting proximity effect and 
ZBA in Andreev conductance are limited by dephasing 
of Cooper pairs due to electron-electron interactions in 
the normal metal. Note that previously Coulomb effects 
in subgap electron transport across NS interfaces were 
studied in a number of workaii"— , however decoherence 
effect of Coulomb interaction was not yet addressed in a 
proper and complete manner. Below we will argue that 
fluctuating electromagnetic fleld produced by fluctuating 
electrons in a disordered normal metal destroys macro- 
scopic coherence of electrons penetrating from a super- 
conductor at a typical length scale L^. The existence 



of this length scale imposes fundamental limitations on 
the proximity effect in NS hybrids at low temperatures 
T < D/L"^. In this temperature range the penetration 
depth of superconducting correlations into the normal 
metal is not anymore given by the thermal length Lt, 
but is limited by the dephasing length which - in con- 
trast to Lt - does not grow at T — 0. We will evaluate 
Andreev conductance G for NS structures in the presence 
of electron-electron interactions and demonstrate that in 
the low temperature limit G essentially depends on L^. 
This dependence allows to directly measure the dephas- 
ing length L^ in transport experiments with NS hybrids. 

It is also interesting to point out that the dephasing 
length Lip derived here for NS systems up to a numer- 
ical prefactor coincides with zero temperature decoher- 
ence length obtained within totally different theoretical 
frameworkiSriS for a different physical quantity - the so- 
called weak localization (WL) correction to the normal 
metal conductance. This agreement demonstrates funda- 
mental nature of low temperature dephasing by electron- 
electron interactions which universally occurs in differ- 
ent types of disordered conductors, including normal- 
superconducting hybrids. On the other hand, as it will 
be explained further below, dephasing of Cooper pairs by 
electron-electron interactions is in several important as- 
pects different from that for single electrons in a normal 
metal encountered, e.g., in the WL problem. 




FIG. 1: (Color online) Hybrid SN structure under considera- 
tion and the diagram describing conversion of a Cooper pair 
into a pair of electrons propagating inside the N-metal. 

The model and formalism. Below we will analyze a 
hybrid SN structure which consists of a normal metallic 
wire of cross-section and length L ^ a attached to 
bulk superconducting and normal electrodes, as shown in 
Fig. f. The contact between the wire and the S-electrode 
is achieved via a small tunnel barrier with cross-section F 
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and resistance i?/ strongly exceeding the wire resistance 
Rj ^ R = L/{aa?)^ where a = 2e^vD is the wire Drude 
conductivity, e is the electron charge and v is the density 
of states per spin direction. 

In order to proceed we will employ the Keldysh 
version of the nonlinear g-niode l^°i^^ adapted to SN 
structures'^^. The effective action for our system defined 
on the Keldysh contour with forward (F) and backward 
(B) parts consists of two terms S = Sw + Sr describing 
respectively diffusive motion of electrons in the wire, 

5„[g,A,$] = '-^i:i[D{dQf ~ AtdtQ + Ai^Q], (1) 
and tunneling between the wire and the superconductor. 



5r[g,A,$] 



Ae'^RiT 



Trr[4c,g], (2) 



where Qsc and Q are taken at superconducting and nor- 
mal sides of the insulating barrier, [x, y\ denotes the com- 
mutator and "Tr" implies the trace over the matrix in- 
dices as well as the integration over times and coordi- 
nates. The covariant derivative is defined as 



<7z 

a. 



(3) 



Here and below Ux.y.z denotes the set of Pauli matrices. 
Both parts of the action ([IJ and ^ depend on the 4 x 
4 dynamical matrix field Q satisfying the normalization 
condition 15{t — t') as well as on the fluctuating 

scalar and vector potentials $(r, t) and A(r, t) which are 
defined on the Keldysh contour and which account for the 
effect of electron-electron interactions. We define = 
± and A± = ^(A^ ± A^) and introduce 
the matrices 
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A+I A-f 
A-f A+i 
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Perturbation theory and Gaussian integration. In what 
follows we will restrict our consideration to energies well 
below the superconducting gap and set 



(T, 



. 

(7, 



5{t-t'). 



(5) 



We will employ the so-called /C-gauge tric k^°i^' which 
amounts to performing the gauge transformation 
g(r, i, t') ^ e'^'£('"'*)g(r, t, t')e-'SK;(r,t') -^^ ^^^^^ ^^an- 
inate linear terms in both electromagnetic potentials and 
deviations from the N-metal saddle point 



U{t^t') 








' sinh(7rT(t-t')) 

-(5(i-i' + 0)i 



I 



(6) 
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This goal is accomplished with the choice of the /C-field 
obeying the following equations 

$+(r,i) = Dd,Kl{v,t) (8) 

-2iDT J dt' coth(7rT(i- i'))5rA^(i",i'), 
$^(r,i) = -Dd,A^{v,t) (9) 

with $K;(r,i) = $(r,0 - dtlC{v,t) and AK.{v,t) = 
A{r,t) — drIC{r,t). After this transformation the ac- 
tion retains its initial form if one substitutes Qsc{t, t') — 

g-<HK;(r,t)g^^(^^^/-)giHK(r,t')^ $ ^ and A ^ A^. 

Treating the tunneling term 1^ perturbatively and 
performing the integration over the g-field, similarly to^^ 
we arrive at the Andreev contribution to our action 



(Trr[gsc,g] Trr [g,e,g])Q. (10) 



The dependence of this term on the electromagnetic 
potentials is encoded both in Qsc and in the average 
of the g-fields. Evaluating Sa within the Gaussian 
approximation we will employ the parametrization22i2i 

Q ^ Qo + iQaoU oW oU - ^QooU oW oW oU with 



W 



ci{r,t,t') di{r,t,t') 

ci{r,t',t) d2{v,t,t') 

di{v,t',t) C2(r,t,t') 

d2{Y,t',t) C2{Y,t',t) 



Here di and ci are respectively the diffuson and the 
Cooperon fields. Expanding the action up to the 
second order in these fields one recovers four different 



contributions S,„ — S-, 



(0,2) 



(1:2) 



(2,1) 



- Sw''^\ where 



5'(',j) contains i-th power of the electromagnetic poten- 
tials and j-th power of W. By direct calculation one 
can verify that the term S'^^'^^ depends only on the diffu- 
son fields which are irrelevant for the problem considered 
here. Hence, our action does not contain the first power 
of the Cooperon fields, and the corresponding propagator 
- the Cooperon C - can be obtained as a solution of a 
linear inhomogeneous equation containing the first and 
the second powers of electromagnetic potentials. 

At this stage we would like to remark that the spin 
structure of the Cooperon analyzed here differs from that 
of the Cooperon encountered, e.g., in the weak localiza- 
tion (WL) problem in disordered normal metals. In- 
deed, representing the Cooperon as a sum of impurity 
ladder diagrams involving retarded (G^) and advanced 
{G^) Green functions one observes that the spin struc- 
ture of the Cooperon responsible for the WL correction 
to the N-mctal conductance is either (ft) or imply- 
ing that both G^ and G^ correspond to either spin up 
or spin down states. In contrast, the spin structure of 
the Cooperon relevant for the proximity induced super- 
conductivity is either (fi) or (If) simply because Cooper 
pairs are spin-singlets. It is straightforward to verify that 
only the Cooperon formed by antisymmetric combination 
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(tJ. — 4,t)/V2 contributes to the subgap Andreev con- 
ductance of NS structures. In the presence of electron- 
electron interactions this Cooperon differs from that cor- 
responding to other possible spin configurations already 
on the level of the first order perturbation theory. 




FIG. 2: First order interaction corrections for the Cooperon. 
The dashed lines represent scattering on impurities and the 
wavy line accounts for electron-electron interactions. 

For illustration let us consider the first order diagrams 
depicted in Fig. 2. While in each of the cases (tt)j (ii) 
and (ti + it)/^/^ these diagrams cancel each other ex- 
actly at T = and in the limit of zero frequencies and 
momentum, no such cancellation occurs for the antisym- 
metric combination (^4- ^ it) I V2 because of extra minus 



sign encountered in this case. Thus, for the latter spin 
combination (which is only relevant here) non-vanishing 
zero temperature dephasing is observed already within 
the first order perturbation theory in the interaction. 

Andreev conductance. Below we will proceed non- 
perturbatively and evaluate the subgap Andreev current 
/ to all orders in the interaction. Defining this current 
as 

I=^J^d^r{SSA/SJC-{r))^ (11) 

and calculating Sa along the lines with the analysisi^ 
which now includes and A^, from Eq. ([TT|) we obtain 



2zye3(i?/r)2 7r 



d^rd^r' / drlra 



{rir,r',T;t)e"'^^) 



sinh(7rrT) 



(12) 

where V = C(r, r; r', 0; t - r/2)e-2»'C+(r.t)+2^K;+(r^^-r/2) 
with the Cooperon C obeying the equation 



{2dr - i$^(r, T - r/2) + i$+(r, T + r/2) - D{dr + iA+{r, T - t/2) + ^A+(r, T + T/2) f) C(r, r; r', r'; T) 

= 6{r~r')S{T-T') (13) 

I 

and V is an external voltage bias. Note that here and be- JC~ which are irrelevant for dephasing of Cooper paira^^. 
low we keep only the fields and /C"*" neglecting $~ and Resolving Eq. we get 



P(r,r',T;i) = 



.^g4K;+(r,i-r)-j/C+(r,t) 



Vxe 



K*+(x(t'),t-(t'+r)/2)-<I.+ (x(i'),t+(t'-r)/2)) _ 



(14) 



x(0)=r' 



What remains is to perform a straightforward Gaussian 
average over ^"'"-fields as well as an average over diffusive 
trajectories. The latter average is handled approximately 
with the aid of the formula (e^)dift — e^^'^'^'^' . As a result 
we find 



(7'(r,r',T;i))$ =2?(r,r';T)e-/('-''-''^) 



(15) 



with /(r, r', r) = Mr, r) + /b(r, r', r) + /t,(r, r', r), 



Mr, r) = - (V++(r, r, 0) - V++(r, r, r)) , (16) 



T i 

fb{r,r',T) J dt J dt' J d'^xdV 



X (V++(x, x', {t - t')/2) - V++(x, x', {t + t')/2)) 
V{r, x; T - t)V{x, x'; t - t')X'(x', r'; t') 
2?(r,r';r) 



, (17) 



where I?(x, x', t) is the diffusive propagator, (r, r' ,t— 
t') = -2i($+(r,t)$+(r',t'))*, and V^^^, V^^ is defined 
analogously. The function /tf, is expressed via the corre- 
lator (/C+$+). We chose to omit it here since ftb remains 
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much smaller than both ft and /;,. 

Eqs. ([n]), (Il5])-([ni) define the central result of our 
work which describes the effect of electron-electron inter- 
actions on the subgap current in diffusive NS structures. 

Quasi- Id structures. Below we will concentrate on 
quasi-ld A^-metal wires (Fig. 1) and set T — a^. In this 
case the differential Andreev conductance G{V) = dl/dV 
takes the form 



G 



^^2 ^(0,0;r) cos(eFr) ^, 
sinh(7rTr) 



/(0.0,r) 



(18) 



with X'(0,0;t) = z?2(0,e-^/^°)/(2La2), where i?2 is the 
second Jacobi theta-function and td — 2L'^ /(tt'^D) is the 
Thouless time. The function / accounts for dephasing of 
Cooper pairs. For ttTt < 1 Eqs. ^ yield 



/(0,0,r)c.-lnf— )+- 

g \TRCJ 



V 



(19) 



In the first term in Eq. p9)) we defined dimensionless 
conductance g = AttvDo^ / L ^ 1 and t^c — RC, where 
C is an effective capacitance. This term is caused by 
spatially uniform fiuctuations of the scalar potential and 
matches with the results^iH. The remaining terms in 
Eq. (jl9p originate from non-uniform in space fluctuations 
in the N-metal and define the new scales in our problem 
- Cooper pair decoherence time = 2'Kva?^/2DT'c and 
decoherence length — yJTJr^^ where Tc ^ 1/vf sets 
a short time cutoff^^"^'' and also T^p ^ tuc- Note that 
up to an unimportant prefactor of order one coincides 
with zero temperature electron decoherence time evalu- 
ated, e.g., for the WL probleroi^"— . 

At this point we would like to emphasize that 
the agreement between the low temperature dephas- 
ing length scales L,^ found here for Cooper pairs and 
previouslyi^^— for single electrons is by no means a pure 
coincidence. Rather this agreement reflects fundamental 
and universal nature of low temperature quantum deco- 
herence caused by electron-electron interactions in vari- 
ous types disordered conductors. At the same time, the 
Cooperon encountered in the WL problem is in many re- 
spects different - both qualitatively and quantitatively - 
from that studied here. As we already indicated above, 
the most important difference is that the spin structure 
of our Cooperon (antisymmetric combination of spin- 
singlets) corresponds to that of a Cooper pair and is 
entirely different from that for the Cooperon in the WL 
problem. In addition, the Cooperon describing propagat- 
ing Cooper pairs in the normal metal is naturally bound 
to the NS interface, which is obviously not the case in 
the WL problem. As a result these two Cooperons are 
defined by formally different diagrammatic series and, 
hence, no apriori conclusions could possibly be drawn for 
our present problem from the Cooperon analysis devel- 
oped for single electrons in disordered metals. 

These differences have several important implications. 
For clarity, let us summarize the most important ones 



again: [i) Unlike single electrons in normal metals. 
Cooper pairs in NS structures get dephased already by 
uniform fluctuations of the scalar potential, as described 
by the first term in Eq. (|19p . (m) unlike in case of the 
WL problem, non-vanishing dephasing of Cooper pairs 
at r = occurs already within the first order perturba- 
tion theory in the electron-electron interactions, see Fig. 
2 and the corresponding discussion above, {Hi) already 
at T = the Cooperon studied here decays differently as 
compared to the Cooperon in the WL problem, cf., e.g., 
our Eq. (19) and Eq. (28) inil, and {iv) at not very low T 
the temperature dependent decay time for the Cooperon 
in NS systems is entirely different from that in the WL 
problem^^. 

o-Rja^G(O) 
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FIG. 3: (Color online) G(0) as a function of L for a = 10 
nm, D = 21 cm^/s. For these parameter values one finds 
1/t^ ^ 0.6 K and L,f, ~ 0.2 ^m. 

Turning to concrete results we first consider the low 
voltage limit eV <^ T. At high temperatures T ^ l/''"<p 
the penetration length of Cooper pairs into the N-metal 
is defined by Lt, while L^p is irrelevant and dephasing 
is only due to spatially uniform fluctuations described 
by the first term in Eq. (|T9|) . In this case the resultsi^ 
are reproduced and one finds G(0) cx At low 

temperatures T l/^c^, on the contrary, Lt becomes 
irrelevant and the penetration length of superconducting 
correlations into the N-metal is set by Then for the 
linear subgap conductance we obtain 



G(0) 



1 f iTRC \ 



s/g 



1 L 



( iTRC ) 
^/2i \ -Tv ) 



2LC(2-f )(22 



8/s 



(20) 

where T{x) is Euler gamma- function and C(x) is Riemann 
zeta-function. The dependence of G(0) on L at different 
temperatures is displayed in Fig. |3l At low T it shows 
a pronounced maximum at i ^ which can be conve- 
niently used for experimental analysis of low temperature 
dephasing of Cooper pairs in NS systems. 

The same information can also be extracted from the 
non-linear subgap conductance G(V) which shows the 
ZBA peak at low voltages^"— . At T -> and L > 
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FIG. 4: (Color online) G{V) at T = and different values of 
L. The parameters are the same as in Fig. 3. 



the width of this peak is roughly determined by ^ ^/t^- 
In particular, for L 3> L^p and T = we get 



r 



(-1) 



(21) 

The non-linear subgap conductance G{V) is depicted in 
Fig. 13] at different values of L. 



Finally we note that our analysis also allows to deter- 
mine the subgap conductance for other geometries. E.g., 
in 3d case the decoherence effect from spatially uniform 
fluctuations is negligible^- and at T ^ I/t'ip the dephas- 
ing of Cooper pairs in the N-metal is controlled by the 
second term in Eq. (|19p with r^p 



In conclusion, we have demonstrated that electron- 
electron interactions yield dephasing of Cooper pairs pen- 
etrating from a superconductor into a diffusive normal 
metal. At low T this phenomenon imposes fundamen- 
tal limitations on the proximity effect in NS hybrids re- 
stricting the penetration length of superconducting cor- 
relations into the N-metal to a temperature independent 
value hip. This new length scale can be probed by mea- 
suring the subgap conductance in NS systems. 

We finally note that our results are qualitatively con- 
sistent with experimental observations'^'* demonstrating 
that the low temperature magnetoconductance of NS 
structures is determined by phase coherent electron paths 
with a typical size restricted by the temperature inde- 
pendent dephasing length L^p rather than by the thermal 
length Lt diverging in the low temperature limit. 
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